Identifying vacancy complexes in compound semiconductors with positron 
annihilation spectroscopy: a case study of InN 
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We present a comprehensive study of vacancy and vacancy-impurity complexes in InN combining 
positron annihilation spectroscopy and ab-initio calculations. Positron densities and annihilation 
characteristics of common vacancy-type defects are calculated using density functional theory and 
the feasibility of their experimental detection and distinction with positron annihilation methods 
is discussed. The computational results are compared to positron lifetime and conventional as 
well as coincidence Doppler broadening measurements of several representative InN samples. The 
particular dominant vacancy-type positron traps are identified and their characteristic positron 
lifetimes, Doppler ratio curves and lineshape parameters determined. We find that Vi n and their 
complexes with Vn or impurities act as efficient positron traps, inducing distinct changes in the 
annihilation parameters compared to the InN lattice. Neutral or positively charged Vn and pure Vn 
complexes on the other hand do not trap positrons. The predominantly introduced positron trap 
in irradiated InN is identified as the isolated Vi n , while in as-grown InN layers Vi n do not occur 
isolated but complexed with one or more Vn- The number of Vn per V"i n in these complexes is found 
to increase from the near surface region towards the layer-substrate interface. 

PACS numbers: 61.72.jd, 78.70.Bj, 71.60.+Z, 61.72.Yx, 68. 55. In 
Keywords: InN, vacancies, density functional theory, positron annihilation 



I. INTRODUCTION 



InN is a significantly cation-anion mismatched semi- 
conductor compound^ with many interesting properties 
and promising applications in opto-, and high-speed elec- 
tronics 2 . Intrinsic point defects have been accounted 
for multiple important mechanisms in the material 3 - - — . 
Amongst them, In and N vacancies and their complexes 
are expected to be the dominant intrinsic acceptors and 
donors, respectively, according to latest density func- 
tional theory calculations^—. Nevertheless, unambigu- 
ous experimental evidence on their nature and charac- 
teristics is still relatively scarce. This is due to limita- 
tions stemming from intrinsic properties of InN, as well 
as challenges related to the growth of high-quality mate- 
rial. Strong surface electron accumulation in polar sam- 
ples^ complicates the fabrication of Schottky contacts 
and therefore the application of standard electrical char- 
acterization methods such as deep level transient spec- 
troscopy (DLTS) and capacitance voltage (C-V) profil- 
ing. Additionally, the high conductivities common for 
early as-grown InN layers together with the unavailabil- 
ity of bulk material strongly limit the use of electron 
paramagnetic resonance (EPR) based techniques. 
Positron annihilation spectroscopy is a powerful method 
for the investigation of vacancy type defects in semicon- 
ductors^ and largely not affected by the above men- 
tioned challenges. Positrons can get trapped and annihi- 
late at neutral and negatively charged open volume sites 
in the crystal lattice due to a locally reduced Coulomb 
repulsion. This increases the positron lifetime and nar- 



rows the momentum distribution of annihilating electron- 
positron (e-p) pairs, both of which can be measured by 
recording the emitted annihilation 7 radiation. 
While direct experimental evidence on the behavior of 
N vacancies has been very limited, previous positron 
annihilation results show that In vacancy (Vi n ) related 
defects are incorporated in concentrations of ~10 16 — 
10 17 cm~ 3 during growth of both molecular beam epitaxy 
(MBE)iSri2, and metal organic chemical vapor deposi- 
tion (MOCVD) 18 InN. Although V In related defect con- 
centrations are found to be low in the 10 16 cm -3 range 
in thick layers of high quality as-grown material^ 7 - the 
experimentally observed concentrations are still by or- 
ders of magnitude higher than what could be expected 
based on first-principles calculations 8 *^. A strong influ- 
ence of the layer thickness on the vacancy concentrations 
has been foundi*, together with a commonly observed 
increase and qualitative change of the vacancy signal 
when approaching the layer-substrate interfac o 7 ' 15 i 18 ' 19 . 
Growth parameters such as polarity* and stoichiome- 
tr y 17 i 18 seem to have only minor impact. This suggests 
that the vacancy formation in InN is not dominated by 
thermal equilibrium processes but rather controlled by 
mechanisms such as local strain, the vicinity to other 
point and extended defects and/or limited surface diffu- 
sion during growth—. Nevertheless, no direct correlation 
between dislocation densities and vacancy formation has 
been found so fa r 17 i 20 . In good agreement with the be- 
havior of negatively charged defects below the branch 
point energy*^, n-type doping of InN by either S i 19 > 22 
or high energy particle irradiation^ leads to an increas- 
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ing incorporation of V\ u related defects, while Mg-doped 
samples (with lowered Fermi level positions) show only 
low concentrations 2 ^. 

Unfortunately, the exact chemical identity of the vacancy 
defects detected in previous studies has remained largely 
unknown. Recent attempts however have shown that 
these limitations could in principle be overcome by care- 
ful modeling of positron annihilation parameters using 
ab-initio calculations 2 ^—. Despite this, a comprehen- 
sive theoretical study of positron annihilation in InN has 
still been missing. 

In the following, we present an extensive identification 
of common vacancy and vacancy-impurity complexes in 
InN combining positron annihilation spectroscopy and 
ab-initio calculations. The employed computational and 
experimental methods are presented in section [ill Cal- 
culated positron lifetimes and momentum distributions 
of annihilating e-p pairs are shown in section IIIII for a 
variety of vacancies and vacancy complexes in InN. In 
section HVl these are compared to experimental data from 
positron lifetime and Doppler broadening measurements, 
and the dominant vacancy defects in positron annihila- 
tion measurements of different as-grown, doped and irra- 
diated InN samples are identified. A critical discussion 
of the results is presented in section IIVDI and possible 
implications are outlined. 



II. METHODS 
A. Experimental 

Positron annihilation measurements are performed us- 
ing a variable-energy (0.5-38 keV) slow positron beam. 
The Doppler broadening of the e-p annihilation radi- 
ation is recorded with two Ge-detectors with a com- 
bined Gaussian resolution function of 1.24 keV (0.66 a.u.) 
and 0.995 MeV (0.53 a.u.) full-width-at-half-maximum 
(FWHM) at 0.511 MeV in the conventional and coinci- 
dence Doppler setup 26 i 27 , respectively. In the latter con- 
figuration, both annihilation photons are detected simul- 
taneously and only counted if energy conservation is ful- 
filled (E tot =1.022MeV). This significantly improves the 
peak-to-background ratio up to 10 6 and sharpens the de- 
tector resolution. To assure statistical reliability, spectra 
of -lxlO 6 and 3xl0 7 counts are accumulated for each 
measurement point in the conventional and coincidence 
setup, respectively. 

In conventional Doppler mode the annihilation peak is 
analyzed using the common integrated lineshape pa- 
rameters 14 which represent the annihilation fractions 
in the low (S) and high (W) momentum parts of the 
spectrum. Integration windows of \pl{S)\ < 0.4 a.u. 
(AE 7 < 0.75 keV) and 1.5 a.u. < \p L (W)\ < 3.9 a.u. 
(2.9 keV < AE 7 < 7.3 keV) are chosen for the S and W 
parameter, respectively. In positron annihilation exper- 
iments the time-integrated annihilation parameter P exp 
(e.g. average positron lifetime, annihilation lineshape, 



S and W parameter) constitutes a weighted sum of the 
characteristic values of the present positron traps Pi and 
the crystal lattice P&, 

Pexp = VbPb + Z?= i ViPi, (1) 

with rjb and rji being the positron annihilation fractions 
of the lattice and the z-th defect, respectively. 
The characteristic defect parameter Pi of a certain 
positron trap can hence be extrapolated from the 
recorded experimental signal P exp , if bulk parameter and 
the positron annihilation fractions are known. The lat- 
ter can be determined by decomposition of the positron 
lifetime spectrum into its n exponential decay compo- 
nents 1 ^, or by a complementary measurement of another 
time-integrated annihilation parameter whose set of char- 
acteristic parameters is well known. The bulk parame- 
ter is usually determined by measurement of a reference 
sample in which no trapping to vacancy defects can be 
observed (confirmed by positron lifetime measurements). 



B. Computational 

Our computational scheme 2 ^ is based on the zero- 
positron-density limit of the two component density func- 
tional theory 2 ^. Valence electron densities are calculated 
self-consistently using the local density approximation 
(LDA) and projector augmented-wave method (PAW) 29 
implemented in the VASP code^S. All electronic struc- 
ture calculations are performed using a 96-atom InN 
wurtzite supercell. Ionic positions are relaxed with a 
convergence criterium of 0.01 eV/A for forces. Indium 
4d electrons are treated as valence and an energy cut-off 
of 400 eV is chosen. The Brillouin zone is sampled with 
a 3 3 Monkhorst-Pack k-point mesh. 
After deriving the electron densities in the lattice, the 
positron densities are solved independently in the cal- 
culated Coulomb potential due to electrons and nuclei 
and the e-p correlation potential. The calculation of the 
positron densities is performed in the so-called "conven- 
tional scheme" , i.e. n + — > is used and the e-p correla- 
tion potential is approximated in the zero positron den- 
sity limit. This is only exact for the case of a completely 
delocalized positron in the crystal lattice (bulk) but is 
also justified for finite positron densities localized at a 
defect site when considering the positron and its screen- 
ing cloud of electrons as a neutral quasi-particle which 
does not affect the surrounding average electron density. 
The LDA2£ in the state-dependent scheme^ 7 - is used for 
the description of many-body effects in the calculation 
of positron annihilation rates and the momentum distri- 
butions of annihilating e-p pairs. In this approximation, 
the momentum density p(p) of annihilating e-p pairs is 
expressed as 

p(p) = 7rr e 2 c£>l J dr e -*P> + (r)^-(r)| 2 , (2) 
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where r e is the classical electron radius, c the speed of 
light and ip + (r), ipj( r ) the positron and electron wave- 
functions, respectively. The summation is performed 
over all occupied electron states and jj is the position 
independent electron-state-dependent enhancement fac- 
tor— which takes into account the increase in annihila- 
tions due to the screening electron density around the 
positron. 

In order to be able to compare the calculated 3D mo- 
mentum density to ID experimental spectra we integrate 
the calculated spectra over the wurtzite m-plane. To 
account for the experimental detector resolution the cal- 
culated momentum distributions are additionally convo- 
luted with a Gaussian resolution function of 0.53 a.u. and 

0. 66 a.u. FWHM for the comparison with coincidence 
and conventional Doppler spectra, respectively. The line 
shape parameters S and W are calculated from the spec- 
tra convoluted with 0.66 a.u. FWHM using the above 
momentum windows. The positron lifetime r is deter- 
mined as the inverse of the annihilation rate A, which is 
given as 

A = — = nrlc J drn + (r)n_(r)7(n_(r)), (3) 

with 7(n_(r)) the enhancement factor— in the LDA. 
In our calculations we made two further simplifications, 

1. e. all defects are considered in the neutral charge 
state and forces due to localized positrons are neglected. 
Adding negative (positive) charge is expected to lead 
to an additional inwards (outwards) relaxation^, while 
the Coulomb force of localized positrons is directed out- 
wards 25 . The exact charge state of defect complexes 
is often unknown in experiments, hence is not safe to 
be assumed without further prove. On the other hand, 
only considering the positron-induced forces 3 -^ leads to 
an overestimation of the outwards relaxation. Neverthe- 
less, careful testing showed that the effects of additional 
charge and localized positrons on the annihilation char- 
acteristics can be neglected compared to changes between 
different defects. Also, trends in the calculated spectra 
are valid irrespective of which approximation is used. 

III. COMPUTATIONAL RESULTS 

The optimized lattice constants for wurtzite InN are 
calculated as a — 3.510 A, c = 1.610 A and u = 0.379, in 
good agreement with literatur e - 11 ' 32 . The positron den- 
sity in the InN lattice is fully delocalized (see Fig. QTa)), 
and a positron lifetime of 157 ps is determined. 
The difference to the experimental lifetime valued 5 - of 
~ 180 ps stems from using the LDA enhancement factor. 
The choice of the LDA is motivated by our focus on the 
calculation of momentum distributions 25 . In this approx- 
imation, good agreement with experiments is achieved 25 
for positron lifetime differences, At = Tdefect — i~b u ik, 
rather than absolute values. When analyzing computa- 
tional positron lifetimes, it should be noted that calcu- 



FIG. 1: (Color online) Iso-surfacc plot of the calculated positron 
density (transparent sphere) in the relaxed lattice structure of the 
(a) InN lattice and (b) 4V N , (c) Vjc n , and (d) 2V In defects. Light 
grey (silver) and dark grey (green) balls indicate N and In ions, 
respectively. 

lated values are sensitive to the size of the relaxed open 
volume (i.e. overlap of positron and electron density), 
and small differences in relaxation of a few percent can 
lead to lifetime differences of a few ps. 
The modeled momentum distribution of annihilating e-p 
pairs in the InN lattice is displayed in Fig. [3J The anni- 
hilation spectrum is symmetrical and is hence folded at 
its center at a.u (511 keV). In the low momentum part 
of the spectrum annihilations with In 5p and N 2p va- 
lence electrons dominate while contributions from tighter 
bound core electrons (In Ad and N 2s) become more im- 
portant in the higher momentum part of the spectrum. 

A. Isolated In vacancy 

The relaxed defect geometry of the V\ n is shown in 
Fig- [3(c) together with the calculated positron density 
in the lattice. An outwards relaxation of the neigh- 
boring N ions from their initial positions is observed in 
good agreement with Duan et ali^. The positron den- 
sity shows a strong localization at the defect site which 
means that V\ n act as efficient positron traps in InN. 
The calculated positron lifetime difference to the InN 
lattice is 85 ps, which is in good agreement with the 
~ 80 ps observed in experiment a 15 ' 19 . Results obtained 
with the atomic superposition (ATSUP) method 3 - 3 - show 
for positrons trapped at the V[ n increased annihilations 
with loosely bound valence electrons (In 5p and N 2p) 
while annihilations with tighter bound core electrons 
(most importantly In Ad) decrease significantly. This 
leads to an overall narrowing of the annihilation peak 
compared to the InN lattice (see Fig. [J). 

For a more detailed analysis of the momentum spec- 
trum we examine its so-called ratio curve in Fig. Eta), in 
which the calculated momentum distribution is displayed 
divided by the spectrum for the defect-free InN lattice to 
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FIG. 2: (Color online) Calculated momentum density of 
annihilating e-p pairs in the delocalized state of the InN lattice 
and trapped at the Vi n , respectively. Trends in the spectra are 
indicated with respect to the InN lattice. 



accentuate the defect-induced changes. The ratio curve 
for the V\ n exhibits a distinct line shape with a maximum 
of roughly 1.08 at the peak center region (0 a.u.). For mo- 
menta above 0.6 a.u. the spectrum drops below 1 and an 
articulate shoulder is visible at 1.2 a.u. which has been 
determined by ATSUP calculations to stem from from 
annihilations with N 2p electrons. At around 3.3 a.u. a 
second broad peak appears with an intensity of around 
0.8 relative to the InN lattice. 

The characteristic relative line shape parameters (for the 
conventional detector resolution of 0.66 a.u.) of the V\ n 
momentum distribution spectrum are calculated as 1.057 
and 0.78 for S and W, respectively. Fig. |4] shows the cal- 
culated values plotted in a SW-plot. 



B. In vacancy complexes 

Recent DFT results^ suggest that Vi n would form 
complexes with each other and predict a significant de- 
crease in formation energies when going from a single 
Vin to 2Vi n and 3Vj n complexes. Motivated by these re- 
sults we calculate the momentum densities and positron 
lifetimes of the 2Vi n and 3Vi n in their most favorable con- 
figurations^. In the 2Vi n pair the two vacancies are lo- 
cated out-of-plane (with respect to the c-plane) on next- 
nearest neighbor positions, and for the 3Vi n complex an 
additional vacancy is added in the c-plane sharing the 
same N atom. As expected, both 2V\ n and 3Vi n com- 
plexes show localized positron densities. Although the 
Vin are only situated on next-nearest neighbor sites their 
open volumes are connected, showing e.g. a handlebar- 
like structure for the positron density at the 2Vi n pah 
(see Fig.[Hd)). 

Positron annihilation characteristics of the 2Vi n defect 
are very similar to the isolated Vi n - The calculated 
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FIG. 3: (Color online) Ratio curves of the calculated momentum 
densities of annihilating e-p pairs in selected vacancy complexes 
in InN. All spectra are convoluted with a Gaussian of 0.53 a.u. 
FWHM (except V^ 6e <*■«■, FWHM = 0.66 a.u.) and divided by 
the momentum density spectrum of the InN lattice. 



positron lifetime is equal and the lineshape parameters 
are very close to the isolated case, with At = 85 ps, 
S = 1.055 and W = 0.77 (see Fig. 0J. In the ratio curve 
in Fig.[3Ka), a slight increase of the peak at 3.3 a.u. is vis- 
ible but the absolute intensity at these momentum values 
is already reduced by ~ 10 3 compared to the peak max- 
imum (see Fig. [2]) and discrimination (in experiments) is 
hence considerably complicated. 
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FIG. 4: (Color online) SW-plot of conventional lincshapc 
parameters S and W for the calculated momentum distributions 
of annihilating e-p pairs trapped at different vacancy type defects 

in InN. 



these defects could in principle promote a localization 
of the positron density. Additionally, the Vn complexes 
are assumed to adapt negative charge states for elevated 
Fermi level positions which additionally supports trap- 
ping. 

To assess whether Vn or its complexes could act as 
positron traps in InN, we calculated the electron and 
positron densities for the isolated Vn and Vn complexes 
of up to 4 N vacancies in the neutral charge state. The 
defect structure and positron density for the 4Vn com- 
plex are shown in Fig. [IJb), exemplarily. We find that 
the positron density for all calculated Vn complexes is 
clearly delocalized and no bound positron state exists at 
Vn or Vn complexes. Hence, isolated Vn and pure Vn 
clusters (in the neutral charge state) cannot be detected 
using positron annihilation spectroscopy. For the case of 
negatively charged Vn complexes the formation of weakly 
localized hydrogenic positron states around the defects is 
expected which possess very similar annihilation charac- 
teristics to the defect free lattice but might be detectable 
through their specific temperature behavior—. 



For the 3Vi n complex we calculate an increase in the 
positron lifetime of roughly 10 ps compared to the iso- 
lated Vin (At — 94 ps). Additionally, the ratio curve 
changes significantly with an increased peak maximum 
and a more pronounced drop at 2.1 a.u. This is also visi- 
ble in the change of the lineshape parameters to S = 1.069 
and W = 0.72. Nevertheless, further analysis shows that 
the relative lineshape of the 3Vi„ and V\ n are very simi- 
lar and the spectrum of the V\ n can nearly be reproduced 
from the 3Vi n spectrum by assuming a positron annihi- 
lation fraction of ij w 0.8 (see Eq. [TJ. This is also vis- 
ible in the SW-plot (Fig. |4} in which the characteristic 
point of the 3Vi n falls on an extension of the line con- 
necting the InN lattice point and the V\ n . For measure- 
ment points lying between the V\ n and the InN lattice, 
the 3Vi„ and V\ n can therefore only be distinguished with 
precise knowledge of the positron annihilation fraction. 
This is in practice limited by the achievable accuracy in 
the separation of lifetime components in positron lifetime 
measurements. 



C. N vacancy complexes 

For most of the Fermi level regions, Vn in InN are 
supposed to be donors in the 3+ charge stated (hence 
repelling positrons) and possess a substantially smaller 
open volume compared to the V\ n , due to the large 
size difference between In and N atoms. Both suggests 
strongly that isolated Vn can not act as positron traps 
in InN. Nevertheless, this may not count for larger com- 
plexes of several Vn which possess a positive binding en- 
ergy (compared to isolated Vn) according to recent first- 
principles calculations 9 . The increased open volume in 



D. Mixed In and N vacancy complexes 

Apart from pure Vi n and Vn complexes, a positive 
binding energy is predicted between isolated Vi n and Vn , 
which should promote the formation of divacancies (Vi n - 
Vn) and larger mixed vacancy complexes in the mate- 
rial^. We calculate the positron densities and positron 
annihilation characteristics for the relaxed lattice struc- 
tures of a variety of mixed vacancy complexes in InN, 
namely Vid-Vn, Vi n -2VN, 2Vi n -VN and Vi n -3VN- As ex- 
pected, the increased open volume leads in all structures 
to a clear localization of the positron density at the defect 
site. The calculated positron lifetimes are close to the Vi„ 
lifetime for the Vid-Vn and Vi n -2VN, with At = 86 and 
88 ps, but increase to 99 and 109 ps for the Vi n -3VN and 
2Vi„-Vn, respectively. Fig. [3fb) shows the ratio curves 
of the computed momentum distributions. A systematic 
trend compared to the isolated Vi n is visible in the spectra 
when adding an increasing amount of Vn around a single 
Vi n . A strong increase of the zero momentum maximum 
to over 1.15 for the Vi n -3VN is visible which is related to 
the increase in open volume. At the same time, the inten- 
sity of the shoulder at 1.2 a.u. decreases with increasing 
number of Vn until it entirely disappears for the Vi n -3VN- 
The ratio curve of the 2Vi n -VN is close to the Vi n -VN for 
lower momentum values but starts to deviate at around 
1.4 a.u. with lower intensities at higher momenta, due to 
reduced annihilation with In 4c? electrons. 
The characteristic line shape parameters of the calculated 
complexes develop accordingly and show a clear deviation 
from the InN - V\ n line in Fig. QJ which is already visible 
for the divacancy and intensifies with increasing number 
of Vn- S and W values for the Vid-Vn, 2Vi„-Vn, Vi n -2VN 
and Vi n -3VN are determined as 1.070 and 0.78, 1.082 and 
0.78, 1.073 and 0.75, and 1.115 and 0.76, respectively. 
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E. Vacancy-impurity complexes 

Ill-vacancy donor-impurity complexes are well estab- 
lished in III- Nitrides. In GaN, a binding energy of 1.6 eV 
has been determined 3 ^ between On and Vg&- The bind- 
ing energy between between Vc a and SiGa is considerably 
smaller (Vc a and Siq & are only next-nearest neighours) 
but still positive 3 ^ Recent ab-initio calculations^ pre- 
dict also in InN a reduction of the defect formation energy 
for Vi n -ON complexes compared to the isolated case, and 
even stronger for the case of Vi n -30^ ■ 
Hence, we take a look at the relaxed defect structures and 
positron densities for the Vi n -ON, V[ n -30N and Fin-iSiin. 
For all complexes the positron density is localized at the 
defect site. The calculated positron lifetimes to the InN 
lattice are again very close to the isolated Vi n case with 
85, 86 and 88 ps for the Vi n -Sii n , %-On and the Vi n - 
30n, respectively. In the ratio curves of the momentum 
densities shown in Fig. [3jc) the peak maximum decreases 
with increasing number of O ions, while the intensity in 
the spectral range above 0.9 a.u. increases including the 
shoulder at 1.2 a.u. and the peak at 3.4 a.u. The form of 
the V[ n -ON ratio curve resembles the case of Vi n trapping 
with a reduced annihilation fraction of r\ f=a 0.8. This can 
also be seen in the SW-plot (Fig. |4]) where the character- 
istic points for all calculated vacancy impurity complexes 
roughly fall on the InN - V\ n line, meaning that their 
lineshape parameters could be reproduced from the Vi n 
values assuming incomplete annihilation fractions. The 
spectrum of the Vi n -Sii n is very close the V\ n and hence 
hardly distinguishable in experiments. The case is dif- 
ferent for the ratio curve of the V[„-30n which possesses 
distinct features with the shoulders at 1.2 and 3.6 a.u., 
respectively, which should be measurable in coincidence 
Doppler measurements. The corresponding lineshape S 
(W) parameters of the calculated vacancy-impurity com- 
plexes are 1.045 (0.83), 1.014 (0.92) and 1.049 (0.77) for 
the Vjn-O-N, V[ n -30N and Vi n -Sij n , respectively. 



F. Summary and discussion 

Our calculations can be summarized with following 
conclusions. Isolated Vn and pure Vn complexes in 
the neutral or positive charge state do not localize the 
positron density and hence cannot be detected with 
positron annihilation spectroscopy. When negatively 
charged, they might act as shallow trapsi^ for positrons. 
Isolated Vi n are efficient positron traps in InN with an- 
nihilation characteristics (positron lifetime, ratio curve, 
lineshape parameters) which are clearly distinguishable 
from the InN lattice values. This is the case for all cal- 
culated defect complexes which include at least one V\ n . 
Apart from the 2Vi n -V^ and V[ n -3VN complexes, their 
calculated positron lifetime values are very similar with 
lifetime differences from 85 ps to 95 ps compared to the 
InN lattice value. Due to the experimental accuracy their 
distinction based on characteristic defect lifetimes is in 



practice not feasible. The case is different for 2Vi n -VN 
and V[ n -3VN complexes, where lifetimes are higher with 
109 and 99 ps difference to the InN lattice, respectively. 
Distinction between the isolated V\ n and the calculated 
complexes complexes can be achieved when considering 
their momentum distributions. For the characteristic 
lineshape parameters of mixed V\ n -nV^ complexes a clear 
deviation from the InN - V\ n line is found in Fig. [¥J with 
a strong increase of the S parameter with nearly constant 
W parameter for an increasing number of Vn around 
the Vi n . Its origin is revealed in the corresponding ratio 
curves as a growing peak maximum and nearly constant 
area for momenta between 1.5 and 3 a.u., where the dom- 
inant contributions to the S respectively W parameter 
stem from. The most distinct feature of these Vi n -nVN 
complexes is identified in the disappearing shoulder at 
1.2 a.u. 

The characteristic points of the calculated vacancy- 
impurity and pure Vj n complexes do not show such devia- 
tion but fall on (an extension of) the InN - V\ n line. Accu- 
rate determination of the annihilation fractions would be 
necessary for an unambiguous identification of the 3V[ n 
complex. The V[ n -30N complex exhibits some unique 
features in its ratio curve, but Vi n -Sii n and 2V\ n are 
hardly distinguishable from the isolated Vj n . 



IV. EXPERIMENTAL SPECTRA AND DEFECT 
IDENTIFICATION 

In order to identify the dominant positron traps in 
common InN material we investigate a variety of as- 
grown and irradiated layers that were grown by different 
growth methods. An overview of a representative selec- 
tion of measured samples is given in Table |U Sample I 
is MBE-grown material 6 which has been irradiated with 
2 MeV He ions to a fluence of 8.9 x 10 15 cm -2 , sample 
II and III are as-grown Si-doped^ and undoped^ 7 . InN 
layers deposited by MBE and MOCVD, respectively. De- 
tails on the growth and characteristics of the samples can 
be found elsewher o 6 ' 7 ' 18 ' 19 ' 36 ' 37 . Sample I and II show 
a strong change in the Doppler broadening signal when 
approaching the interface, which is a common feature 
in several previously investigated InN sample s 7 ' 15 ' 18 ' 19 . 



TABLE I: Studied samples and experimentally determined 
lineshape parameters of the respective dominant positron trap. 
All parameters are extrapolated to saturation trapping and 
divided by the reference values measured in InN lattice. 



Layer: Interface: 



ID 


Sample 


S/S re f 


w/w ref 


S/S r ef 


w/w rcf 


I 


MBE, irr. 


1.042 


0.80 


1.083 


0.78 


II 


MBE, Si-doped 


1.051 


0.83 


1.077 


0.81 


III 


MOCVD 


1.052 


0.81 


/ 


/ 
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0.70 



1.00 



1.02 



1.04 



1.06 



1.08 



1.10 



s/s, 



FIG. 5: (Color online) Experimentally determined S and W 
values from conventional Doppler broadening setup for different 

samples (see Table [Xp . Open symbols correspond to directly 
measured values, closed symbols have been extrapolated using 
trapping fractions estimated from positron lifetime measurements. 

Light colors are used for points at high implantation energies, 
close to the interface (see arrow), while dark colors indicate points 
for low implantation energies, in the InN layer region. 



Therefore, the interface region is investigated separately 
in these two samples (see section H"VC[) . 
Conventional Doppler broadening (Fig. [S]) as well as co- 
incidence Doppler (Fig. [6]) spectra are recorded for all 
samples. The recorded momentum distributions for sam- 
ples I and II are extrapolated to saturation trapping, 
using annihilation fractions determined from previously 
recorded positron lifetime measurements. An extensive 
description of the positron lifetime data is published else- 
where^^. All experimental momentum distributions are 
shown as ratio curves divided by the spectrum of a suit- 
able reference for the InN lattice. The reference sample 
has been carefully analyzed using positron lifetime spec- 
troscopy and no positron annihilation in trapped states 
at vacancy defects is observed^. All samples are mea- 
sured perpendicular to the c-axis. 



A. Irradiation-induced defects 

Positron lifetime measurements^ of the layer region of 
sample I show one dominant positron trap with a charac- 
teristic lifetime of 260 ps, and an annihilation fraction of 
rj = 0.47 is determined for the trap. The recorded line- 
shape parameters from conventional Doppler broaden- 
ing measurements are extrapolated accordingly and dis- 
played in Fig. [3] Characteristic line shape parameters of 
S = 1.042 and W = 0.80 are evaluated for the dominant 
defect in this region. 

The extrapolated ratio curve of sample I is shown in 
Fig. Ufa). When comparing the experimentally deter- 
mined defect ratio curve to the calculated momentum dis- 



tributions in section lll Bl we find good agreement with the 
spectrum of the isolated Vi n for the most of the spectral 
range. In the central region of the peak slightly higher in- 
tensities are found in the calculated spectrum compared 
to the experimental one. This region is mostly sensitive 
to the size of the open volume of the positron trap, with 
higher intensities for larger volumes. The calculated 2V\ n 
and V[ n -ON complexes possess an overall rather similar 
shape of the ratio curve and based on this also have to be 
considered as possible sources of the defect signal. Nev- 
ertheless, the available lifetime data does not support the 
reduced annihilation fraction which would be required for 
an identification of the experimental spectrum with the 
Vui-On- Although the differences between the Vj n and 
2Vi n complex are more subtle and hence the 2Vi n cannot 
be ruled out contributing to the signal, the ratio curve 
of the Vi n gives the best overall approximation of the ex- 
perimental spectrum. Therefore we identify the positron 
trap created in irradiated InN with the isolated V\ n . 



B. Defects in as-grown samples 

The experimentally determined S and W parameters 
from conventional Doppler broadening for the as-grown 
samples II and III are displayed in Fig. [SJ Positron 
lifetime measurements^ show a annihilation fraction of 
r/ = 0.34 in the dominant positron trap for sample II. 
Based on this, the characteristic lineshape parameters 
for this defect are estimated as 1.051 and 0.83 for S 
and W, respectively. This is close to the as-measured 
line shape parameters of sample III, i.e. S = 1.052 and 
W = 0.81. For this sample no lifetime data is avail- 
able. In Fig. HJb) the respective ratio curves recorded 
with coincidence Doppler measurements are displayed. 
The extrapolated ratio curve of both sample II and the 
as-measured spectrum of sample III show a very similar 
line shape. The bigger scatter in the former is due to a 
smaller annihilation fraction. 

Compared to sample I, the as-grown samples II and III 
show several differences in their ratio curves. First, the 
intensity in the peak center region is clearly increased. 
The intensity difference to the InN lattice is thereby mag- 
nified by about 35 % compared to the spectrum of sam- 
ple I. This is supported by very accurate statistics in 
this spectral region. Second, a significant decrease of the 
shoulder at 1.2 a.u. is visible, also with high statistical 
accuracy. Third, the drop at 2 a.u. is less pronounced, 
followed by slightly higher intensities in the high mo- 
mentum region of the spectrum. Nevertheless, stronger 
scatter starts to dominate this region. 
A comparison with the calculated defect spectra in sec- 
tion IIIII reveals that these changes coincide with the ef- 
fects of the decoration of a Vj n by Vn, as presented in 
Fig. [3fb). Especially the characteristic decrease of the 
shoulder at 1.2 a.u. in ratio curves of the experimental 
spectra cannot be correlated with any other calculated 
vacancy defect complex (see additionally Hautakangas 
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^ Sample III (layer) 
(y Sample II (layer) 
O Sample I (layer) 



V -V 

In N 




Sample I (interface) 
\/ Sample II (interface) 
O Sample I (layer) 



V. -3V. 




(c) 



12 3 4 

Momentum (a.u.) 

FIG. 6: (Color online) Experimental coincidence Doppler 
spectra of the investigated samples in the layer (a,b) and interface 
(c) region. The data has been divided by a suitable reference 
spectrum for the InN lattice. Computational ratio curves are 
shown for comparison. 



et alM). This is also expressed in the observed devia- 
tion of the characteristic lineshape parameters of sample 
II and III from the line determined by the characteris- 
tic points of sample I and the InN lattice in Fig. [5] A 
similar trend can be observed for the calculated V\ n -nV^ 
parameters in Fig. 2) Judging from the amount of the ob- 



served changes in the ratio plots an identification of the 
experimental spectra with V[ n -VN is most feasible, with 
possible influence from the V\ n -2V-^. This assignment is 
in good agreement with the positron lifetime data 19 . 



C. Defects at the Interface 

A strong change in the Doppler broadening signal is 
observed for sample I and II close to the interface region, 
as visible in Fig. [5] From positron lifetime data in this re- 
gionLiii we determine a positron annihilation fraction of 
r\ = 0.47 and rj — 0.52 at the dominant positron traps in 
samples I and II, respectively. The characteristic defect 
lifetime at the interface is slightly increased compared to 
the layer region. The extrapolated line shape parame- 
ters are determined as 1.083 (0.78) and 1.077 (0.81) for S 
(W) in samples I and II, respectively, and show an even 
stronger deviation from the line defined by the character- 
istic points of sample I and the InN reference as already 
observed in section HV Bl The characteristic SW points 
of both interfaces coincide within the statistical accu- 
racy. This indicates that the same dominant positron 
trap is present in both samples, which is also shown in 
the extrapolated ratio curves shown in Fig. E^c). In both 
samples a strong increase in the peak center intensity to 
about 1.12 is visible, which is over 2-fold compared to 
that observed in the irradiated layer. Additionally, the 
signal drops straight to the minimum at 2 a.u. without 
showing anymore the shoulder which is visible in the layer 
region of both samples. The observed trends are qualita- 
tively very similar to the ones described in the previous 
section for the layer region of samples II and III, but 
intensified. Therefore, we identify the induced changes 
with an increase in the decoration of Vi n with Vn ■ When 
comparing to the calculated momentum distributions in 
section IIIH best agreement is found for the spectrum of 
the Vi n -3VN complex. Results from both conventional 
Doppler broadening (Fig. [5]) and positron lifetime spec- 
troscopy additionally support this assignmen t 19 ' 38 . 



D. Summary and discussion 

Based on the above presented data we are able to iden- 
tify the dominant positron trap created in high-energy 
particle irradiation of MBE grown InN layers as the iso- 
lated Vi n , while in as-grown MBE and MOCVD mate- 
rial the observed defect is a mixed V[ n -VN vacancy com- 
plex. The changes at the interface in both irradiated 
material and as-grown layers are assigned to the for- 
mation of larger Vi n -nV^ complexes with an average of 
about 3 Vn surrounding the V\ n . It has to be noted that 
the positron implantation profile for the chosen interface 
point (12 keV) is already considerably broadened and 
hence signal from a wide area of the layer is averaged. 
The observed complexing of V\ n and Vn is in good agree- 
memei with existing data on the electron mobility in 
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these samples^ and might be the reason behind observed 
superior transport properties. 

The measured Doppler broadening of the annihilation 7 
radiation is caused by the momentum component of anni- 
hilating e-p pairs along the detection direction^. All our 
experimental spectra arc recorded perpendicular to the 
c-axis but the exact detection direction in the c-plane is 
not identified. We compare the experimental spectra to 
computational momentum distributions integrated over 
the m-plane. The induced uncertainty could in principle 
cause difficulties for the exact defect identification, as 
small differences between the different lattice directions 
are present^. Nevertheless, these are minor compared 
to differences between the calculated spectra of different 
defect complexes and hence pose no problem in our case. 
This counts especially when interpreting changes in the 
momentum distributions rather than absolute values. 
Our data accentuates the advantage of high-quality co- 
incidence Doppler— spectra for the identification of de- 
fect identities. Apart from the advanced experimental 
resolution compared to conventional Doppler measure- 
ments which helps to distinguish defect-induced changes 
(see Fig. [3{a)), important features in the spectra might 
be lost when only regarding the integrated lineshape pa- 
rameters S and W. This is the case, e.g., for the observed 
disappearance of the shoulder at 1.2 a.u. which is identi- 
fied as an unique feature of the Vi n -nVN complexes. The 
commonly used integration windows for the S and W pa- 
rameters cannot include this momentum area, in order 
to avoid correlation effects^. Irrespective of this, a shift 
of the lower limit of the W-parameter window to 2.2 a.u. 
(for a detector resolution of 1.24 keV) would be beneficial 
in future positron experiments in InN in order to maxi- 
mize the visibility of defect-induced changes (see Fig. [3J . 



V. CONCLUSION 

We present the identification of dominant vacancy-type 
positron traps in different representative InN samples 
by combining positron lifetime and Doppler broadening 
spectroscopy with ab-initio calculations of the positron 
annihilation characteristics. Calculated momentum dis- 
tributions of annihilating electron positron (e-p) pairs in 
different vacancy type defects are compared to high reso- 
lution measurements of the Doppler broadened e-p anni- 
hilation radiation acquired with the coincidence Doppler 
technique. We find that isolated Vn or pure Vn complexes 
do not trap positrons. During high-energy particle irra- 
diation isolated V\ n are created, which are the dominant 
positron trap in such samples at room temperature. In 
as-grown InN samples on the other hand, V\ n do not oc- 
cur isolated but complexed with one or more Vn- The 
observed changes of the Doppler broadening signal close 
to the interface are identified with an increasing amount 
of Vn surrounding the Vj n . The characteristic lifetimes, 
ratio curves and S and W parameters for the identified 
defects are determined. 
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